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Abstract 



The neutrinoless double beta decay (O^/3/3-decay) transition to the first ex- 
cited + collective final state is examined for A=76, 82, 100 and 136 nuclei 



by assuming light and heavy Majorana neutrino exchange mechanisms as 
well as the trilinear R-parity violating contributions. Realistic calculations of 
nuclear matrix elements have been performed within the renormalized Quasi- 
particle Random Phase Approximation. Transitions to the first excited two- 
quadrupole phonon + state are described within a boson expansion formalism 
and alternatively by using the operator recoupling method. We present the 
sensitivity parameters to different lepton number violating signals, which can 
be used in planning the Oi^/5/5-decay experiments. The half- life of Oi^/3/3-decay 
to the first excited state Of is by a factor of 10 to 100 larger than that of the 
transition to the ground state, 0+ . 
PACS numbers:23.40.BW,23.40.HC 



I. INTRODUCTION 

The neutrinoless double beta decay (O^/3/3-decay) , which violates the total lepton number 
by two units, is the most sensitive low-energy probe for physics beyond the standard model 
(SM) |'[|-§]. The observation of the Oz//3/3-decay would provide unambiguous evidence that 
at least one of the neutrinos is a Majorana particle with non-zero mass B, This conclusion 
is valid without specifying which from the plethora of possible 0^/?/?-decay mechanisms trig- 
gered by exchange of neutrinos, neutralinos, gluinos, leptoquarks etc is the leading one. The 
current experimental upper limits on the Ou (3 /3-decay half-life impose stringent constraints, 
e.g., on the parameters of Grand Unification (GUT) and supersymmetric (SUSY) extensions 
of the SM. 

There is a continuous, both experimental and theoretical, activity in the field of Qv/3/3- 
decay. An interesting issue is what are the implications of the neutrino oscillation phe- 
nomenology to the O^/3/3-decay. We note that the results of the solar 0, atmospheric || and 
terrestrial || neutrino experiments provide a convincing evidence of neutrino oscillations, 
which require non- vanishing masses for neutrinos as well as neutrino mixing [|IIJ]. The neu- 
trino oscillations are sensitive to the differences of the masses squared and cannot distinguish 
between Dirac and Majorana neutrinos. Nevertheless, if assumptions about the character 
(Dirac or Majorana neutrinos), the phases and the neutrino mixing pattern are considered, 
one can derive estimates for the effective Majorana electron neutrino mass < m u > responsi- 
ble for O^/3/3-decay. The current viable analysis implies the effective neutrino mass < m v > 



to be within the range 10 eV < < m u > < 1 eV [p!l] , P^1 . The present generation of 



0^/9/9-decay experiments |T3|-[15| achieve sensitivities of T®Y 2 ~ 10 24 — 10 25 years for different 



isotopes ||l3Hi5ll. It corresponds to < m u >m 0.5 — 1. eV Hlq , |l7|l , i.e., these Ou/3/3-decay 
experiments allow already to discriminate among various neutrino mixing schemes. 

The neutrino oscillations imply that perhaps we are close to the observation of the 0vf3f3- 
decay. This would be a major achievement. Maybe it is enough to increase the sensitivity 
to < m v > by about an order of magnitude, i.e., the 0z//3/3-decay experiments should be 



sensitive to half-lifes of 10 27 — 10 28 years for the ground to ground transitions. We hope 
that these data will stimulate new experimental activities. Ambitious plans are underway 
to push the upper constraints on lepton number violating parameters further down. By 
using several tons of enriched 76 Ge, the GENIUS experiment is expected to probe < m u > 
up to 10 -2 eV [IS|. The CUORE experiment intends to search for rare events with the 
help of a cryogenic Te02 detector with high energy resolution ||19|| . The ongoing NEMO 3 
experiment, now under construction in the Frejus underground laboratory, will measure up 
to 10 kg of different double beta decay isotopes 1 20f| . Both CUORE and NEMO 3 have a 
chance to reach a sensitivity to the effective neutrino mass < m u > in the order of 0.1 eV 



m 

It is worth to examine also other possibilities to increase the sensitivity of 0z//?/3-decay 
experiments. Till now, the attention was concentrated mostly to the 0z//?/5-decay transition 
to the ground state of the final nucleus. However, there might be a chance that the transitions 
to the excited + and/or 2 + final states are more favorable experimentally, at least for a 
particular mechanism for the 0^/3/3-decay |T^f2T|. Generally speaking, transitions to the 
excited states are suppressed due to the reduced Qpp value. However, this restriction can 
be compensated by a possible lower background due to a coincidence of the (3 particles with 
the 7 or 7's from the excited final state. The possible advantage depends on the ratio of 
corresponding nuclear matrix elements to the excited and to the ground state. If their values 
are comparable, the Oz//3/3-decay experiment measuring transitions to ground and excited 
final states could be of a similar sensitivity. 

Recently, this issue received an increasing attention from many experts in the field. The 
0i//3/3-decay of 76 Ge and 100 Mo to the first excited 2f final state has been investigated in Ref. 
|22|| by assuming massive Majorana neutrinos and right-handed weak currents. It was found 
that the zero neutrino double beta decay transition probabilities for the + — > 2 + decay 
are strongly suppressed due to higher partial waves of the emitted electrons needed in the 
+ — > 0+ s transition. But in the case of the O^/3/3-decay to lowest and first excited + states 
the two emitted electrons are preferentially in a S1/2 wave. Therefore, this decay channel is 
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more favored. Recently, the first realistic calculation for + — > 0^ decay has been performed 
for A=76 and 82 nuclei within a higher Quasiparticle Random Phase Approximation (QRPA) 
||23|| . The 0^/?/9-decay mechanisms mediated by light Majorana neutrinos within the left- 
right symmetric model have been discussed. The conclusion was that the transition to 
excited collective + final states are reduced compared with the decay to the ground state. 
It would be worthwhile to test this result also within other nuclear approaches. 

The aim of this work is to examine the O^/3/3-decay of 7e Ge 82 Se, 100 Mo and 136 Xe to 
the first excited 0^ state. We shall discuss the mechanisms induced by light and heavy 
Majorana neutrino exchange as well as those with trilinear R-parity violation. The nuclear 



matrix elements will be evaluated within the renormalized QRPA (RQRPA) p4j , 25| , which 
take into account the Pauli exclusion principle. We shall also consider the contributions 
to O^/3/3-decay coming from the momentum-dependent induced nucleon currents, which has 
been found to be significant for ground state to ground state transitions [|16|]. We note 
that this contributions were ignored in a similar study 1 123|| . The collective two-quadrupole 



phonon state 0{ , will be described by two different approaches proposed in Refs. []2B| and 
]2"7^P§|, respectively. Finally the sensitivity parameters for a given isotope, associated with 
different lepton number violating signals and nuclear transitions, will be calculated. Also a 
discussion about possible projects of future experimental searches for 0z//3/3-decay to excited 
final states + , will be presented. 

The paper is organized as follows. In section II, basic formulae relevant to the 0vf3f3- 
decay mechanisms, are presented. In section III two approaches meant to calculate the 
transitions to ground and excited states, are described. In section IV we calculate the 0v(3/3- 
decay matrix elements for A=76, 82, 100 and 136 nuclei via RQRPA. We also determine 
the sensitivity of transitions to the collective excited state Of, to lepton number violation 
associated with the Majorana neutrino mass and R-parity breaking. The perspectives of 
measuring + — > Of decays are analyzed. The summary and final conclusions are presented 
in Sec. IV. 



II. THE Oi//3/3-DECAY HALF-LIFE 

The theory of light and heavy Majorana neutrino mass modes of O^/3/3-decay have been 
reviewed, e.g., in Refs. 0J^,|,|T^]. The trilinear R-parity violating mode of 0z//3/?-decay has 
been presented in Refs. pOHSBf . Without going into details of derivations, we summarize 
the basic ingredients of these modes of Oz//3/3-decay. 

A. Majorana neutrino mass mechanism 

The half-life of 0^/3/?-decay associated with light and heavy Majorana neutrino mass 
mechanism is given as 

<f 77? ^> 

[2ft]- 1 = G 01 \ — M l itl> + V N M^\ 2 . (2.1) 

Here, Gqi is the integrated kinematical factor [l,34]]. The lepton number non-conserving 



parameters, i.e., the effective electron Majorana neutrino mass < m v > and 77 , are given 
as follows: 

<m u > = E (L&) a ^ m k , Vn = E (^) 2 2 fc -^, (2.2) 

1 1 * 

with irtp (m e ) being the proton (electron) mass. U L is the unitary mixing matrix connecting 

left-handed neutrino weak eigenstates uil, to mass eigenstates of light Xk and heavy Nk 

Majorana neutrinos with masses m^ (m^ << 1 MeV) and M^ [M^ » 1 GeV), respectively. 

We have 

"il= E UfcxkL+ E U&Nu, (l = e,HT). (2.3) 

k=light k=heavy 

rp 

Vki Nk satisfy the Majorana condition: v^h — C V^, N^k — C N k , where C denotes the 
charge conjugation while £, S are phase factors; the eigenmasses are assumed positive. 

The nuclear matrix elements associated with the exchange of light (M<m„>) an d heavy 
neutrinos (M^ eavy ), including contributions from induced nucleon currents can be written 



as a sum of Fermi, Gamow- Teller and tensor components ||16|| : 



M K = f + M£ T + M% {K =< m v >, Vn ). (2.4) 

9a 

Here, g& = 1.25. Using the second quantization formalism, M.^ can be expressed in terms 
of relative coordinates as follows: 

{) v In " 
jn> jp> J 

x(j ) (l),p'(2);J\f(r 12 )r+r+0 )C (12)f(r 12 )\n(l),n'(2);J) 

x(0+ || {cp n ,]j || rmfHrmfirm^irrm || [c + c n ]j || Of). (2.5) 

Here, f{r 12 ) is the short-range correlation function [f| and Ojc(12) (/C =< m u >, rj N ) 
represents the coordinate and spin dependent part of the two-body Oi//3/3-decay transition 
operator 

O k (12) = - Hp ^ l2) + H% T (r 12 )a 12 + H$(r 12 )S l2 . (2.6) 

9a 

Also, the following notations have been used: 

_ i | o _ r i2 
ri2 — r l — r 2, ?"12 — r 12 , ri2 — — , 

Sl2 = 3(0 ! i -fl2)(<?2 -fl2) -0"12, 0"12 = ^1 ' 0*2- (2-7) 

where ri and r 2 are coordinates of the beta decaying nucleons. The radial part of the light 
and heavy neutrino exchange potentials Hf miJ> (r 12 ) and if / JV (r 12 ) (/ = F, GT, T) can be 
written as 

r <m v> , ^ _ 2 R [°° sm(qr 12 ) 



with 



1 2 R r°° 

H V I N (r 12 ) = J— / smiqn^hriq 2 ) q dq (2.8) 

m p m e 7r^ ri2 ^o 



2, __ 2, 2^ , l^V 2^(g 2 )^(g 2 )g 2 , 2^(g 2 )g 2 



^ ) = <fc(? ) + 3 4m , 3 2mp ■ 3 4m 2 > 

. ,*. ^ 9A{q 2 )9p{q 2 )q 2 1 sfW ! <&(g 2 )g 2 ,„ Q , 

hAq ] ~ 3 2^ S^ 2 "" + 3 4m 2 " (2 ' 9) 



Here, R = r Q A 1 ^ is the mean nuclear radius [lfj with r = 1.1 fm. E l , E? and E m [J) are 



the energies of the initial, final and intermediate nuclear states with angular momentum J, 
respectively. The momentum dependence of the vector, weak magnetism, axial-vector and 



pseudoscalar formfactors (gv(q ), 9m(q ), 9a(q ) and gp(q )) can be found in Ref. |L6 . 

We note that the overlap factor (J n mf\J' 7T m i ) and the one-body transition densities 
< .Prtii || [c+c n ]j || 0+ > and < 0^ || [c+c n ]j || J n rrif > entering Eq. ( |2.5| ) must be 
computed in a nuclear model. 

B. The trilinear R-parity violating mechanisms 

The minimal supersymmetric standard model (MSSM), which is the simplest extension 
of the SM, preserves the R-parity, i.e., also the total lepton number. We remind that im- 
parity is a discrete multiplicative symmetry defined as R p = (— i) 3B + L + 2S, 5 where S,B and 
L are the spin, the baryon and the lepton quantum number. Thus R p — +1 for SM particles 
and Rp — — 1 for superpartners. 

Generally, one can add trilinear R-parity violating terms to the superpotential of MSSM 

W % = \ %]k UL 3 E% + K 3k UQ 3 Dl + vijLjHi + KjkU'Dpt, (2.10) 

where i,j, k denote generation indices. Here L, Q stand for the lepton and quark doublet 
left-handed superfields while E c , U c , D c for the charge conjugated lepton, up and down 
quark singlet superfields, respectively. The terms proportional to A and A' violate the lepton 
number while those proportional to A" violate the baryon number. 

The 0z//3/3-decay can be induced by different trilinear R-parity violating mechanisms, 



partially determined by different products of the parameters A and A' | 12]j33| ) |35 ]. Here, we 
consider those mechanisms which lead to the most stringent constrain on the A'm parameter. 
They are triggered by exchange of gluinos and neutralinos. The corresponding Feynman 



diagrams can be found, e.g., in Ref. [[L7|,[33|]. If the masses of the SUSY particles are assumed 



to be of about the same value, there is a dominance of the gluino-exchange mechanism [33 . 



This conclusion is expected to be valid also for the O^/3/3-decay transitions to excited + 
states. 

The Ozz/S/S-decay half-life, associated with exchange of gluinos, is HJSc 



[T 1/2 (0+ - 0+)]- 1 = Goi h + ^s\ 2 \M x/lll \ 2 . 



(2.11) 



The effective $ p violating parameters r\g and rj^ can be expressed by means of the fundamental 
parameters of the MSSM, as follows: 
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Here, a s = gf/(47r) is the SU(3) C gauge coupling constant. ma L , m d - and m g are masses of 
the u-squark, d-squark and the gluino. 

At the level of hadronization, the dominant mechanism is the pion realization of the 
underlying AL = 2 quark-level Oz//5/5-transition dd —* uu + 2e~ 



The nuclear matrix 



element M\' in can be written as a sum of contributions originating from one and two 
pion-exchange modes. Thus, we have pB 



Mx'm 



M 1 * + M 

2 



2- 






r27T 



(2.13) 



where 



with 



R 



M% T T = (0+\J2T?T+a l -a J Fg>(m n r 7T ) — \0t), with fc = l,2 
¥i rij 



(2.14) 



F«(x) 



e-, 4 1} (x) 



(3 + 3x + x 2 



F$(x 



ar 
(x) = (x - 2)e _!B , i^ 2) (x) = (x + l)e" 



(2.15) 
(2.16) 



Here, m^(= 850 MeV) and m n are the mass scale of nucleoli formfactor and the mass 



of the pion, respectively. Values of the structure coefficients a (k=l,2) are J33|: a 



,1- 



-4.4 x KT 2 , a 2n = 2x 10 _1 . 

Having in mind the forthcoming calculations within the RQRPA, it is useful to rewrite 
Adx'm m ^he f° rrn given by Eqs. ( |2.5| ) and ( |2.6| ) (/C = A'm). One finds that the Fermi part 
of the pion-exchange potential is equal to zero (Hp 111 ^^) = 0) and the the Gamow- Teller 
and tensor parts are given by: 

/ \ 2 /A 



2 



YYv m J TTvg \ o 



^' 111 (r 12 ) = — ^ f^^F^(m w r 12 ) +a 2 -F^(m.r 12 )) . (2.17) 

Vm p / m e V3 / 

III. ONE BODY TRANSITION DENSITIES 

The nuclear matrix element M.^ (/C =< m u >, r/ N A'm), in Eq. ( |2.5| ), is associated 
with the 0z//3/3-decay to the ground state, 0^ s , or to any of the excited + states in the 
final nucleus. Its evaluation requires the description of the initial (|0^)), final (|0j)) and the 
intermediate states (all in different nuclei) with angular momentum and parity J 77 (| J^rriij)), 
within a given nuclear model. Then, the one-body transition density, entering the expression 
for M-ic-i can be calculated and consequently the chosen matrix element is readily obtained. 

The standard QRPA (based on the quasiboson approximation) and the RQRPA have 
been intensively used to calculate nuclear matrix elements for the double beta decay 
P JT6|j23| - p5| , p2| -|34]| . The RQRPA includes anharmonicities (the ground state is less cor- 



related than in the standard QRPA) and there is no collapse of its first solution within 
the physical range of the particle-particle interaction strength. Within schematic models, 
it has been shown that by including the Pauli Exclusion Principle (PEP) in the QRPA, 
good agreement with the exact solution of the many-body problem can be achieved even 
beyond the critical point of the standard QRPA pB]. The RQRPA takes into account the 



PEP in an approximate way. Nevertheless, it is enough to avoid the main drawback of the 

9 



standard QRPA and to reduce the sensitivity of the calculated observables to the details of 
the nuclear model. The RQRPA has been used in our previous studies of the double beta 
decay P JT6| . [29|J3^J33|J37[] . Here we apply this approach to calculate the Oz//3/3-decay to first 



excited states + . 

The final nuclei for A = 76, 82, 100 and 136 double beta decaying systems are 76 Se, 82 Kr, 
100 Ru and 136 Ba, respectively. The first excited Of state of these nuclei is believed to be 
member of the vibrational triplet + , 2 + , 4 + . This state can be described as follows 

|o 1 -) = ^ f {r^ t ®r^}0|o+j, (3.i) 

where T 2 is the creation quadrupole phonon operator. The experimental energies of Of 
(E(0f)) and 2+ (E(2f)) states relative to the ground state energies are 

[E(0f),E(2+)] = [1.122,0.559] MeV for A = 76, 
= [1.488, 0.777] MeV for A = 82, 
= [1.130,0.540] MeV for A = 100, 
= [1.579, 0.818] MeV for A = 136. (3.2) 

One notices that the energy of the 0^ excited state is about twice the energy of the 2f 
excited state. 

The nuclear states of interest are described as charge changing (pn- RQRPA) and charge 
conserving (ppnn-RQRPA) modes of the RQRPA approach. 

In the framework of the pn-RQRPA, the m th excited state of the intermediate odd-odd 
nucleus, with the angular momentum J and the projection M, is created by applying the 
phonon-operator Qj M -n on the vacuum state |0^ PA ): 

\m, JM*) = Q7LKpa) with QImA^rpa) = 0- (3.3) 

Here \0~^ PA ) is the ground state of the initial or the final nucleus and the phonon-operator 
Qj M -k is defined by the ansatz: 

Q%> = E fe,j^ V, JM) + Y^ Jn A(pn, JM)] . (3.4) 

pn 

10 



-^Tpn, •/*-)> Y(vn,j^ denotes free variational amplitudes, which are calculated by solving the 
RQRPA equations. 

The first excited 2 + state of the daughter nucleus is assumed to have one quadrupole- 
phonon character. Within the ppnn-RQRPA (allowing for two proton and two neutron 
quasiparticle excitations only) this state is defined by: 

|2+) = r^ + |0+ PA ) with I1 M+ |0+ PA > = 0, (3.5) 

where 



Am+ = E [ J R( P y, 2+) ^ t (p,p',2M) + ^ p , 2+) i(p,p / ,2M) 



+ 



p<p' 



]T [i?f n;n , j2+) At(n, n', 2M) + S^ 2+) A(n, ri, 2M)] . (3.6) 

n<n' 

A\tt' , JM) and A(tt', JM) (r = p,n and r' = p',n') are the two quasi-particle creation 
and annihilation operators coupled to the good angular momentum J with projection M 
respectively, defined by: 

A (TT ,JM) = j 2 2_^ ( ~'j T m T j T ,m T , a Tm T a T'm' T i 

v * tt ) m T ,m' T 

A(tt', JM) = (a\tt', JM)V . (3.7) 



The vacua defined by Eqs. (3.5) and (3.3) are in principle different from each other. However 
the differences induce corrections to the matrix elements considered, of higher order and 
therefore they are neglected. The quasiparticle creation and annihilation operators (a+ m 
and a TniT , r = p,n) have been defined through the Bogoliubov- Valatin transformation 

u T v T \ /c+ m \ 

(3.8) 
v T u T j y c TrnT J 

where c+ m (c rm ) denotes the particle creation (annihilation) operator acting on a single par- 
ticle level with quantum numbers (n T , l T ,jr)- The parameters u, v are occupation amplitudes 
and the tilde symbol indicates the time-reversal operation, e.g. a TniT = (—l)^ r ~ mT a T _ mT . 
Let us now denote by T> pn and T> TT > (r = p, n) the following expectation values: 

11 



(0+ PA \[A{pn, JM), A\p'n\ JM)}\0+ PA ) = 5 pp ,5 nn/ V pn , 
(0^|[A(TT^JM),^(^^JM)]|0+ PA ) = (5 T Av-(-l) JV+ ^" J ^5 T gP TT ,. (3.9) 

Here, the exact expressions of the commutators are taken into account. The calculation of 
V factors is discussed in Ref. [ p5| , p8[ . 

Solving the pn-RQRPA (ppnn-RQRPA) equations, one gets the renormalized amplitudes 



X, Y (i?, S) with the usual normalization XX — YY = 1 (RR — SS = 1). They are related 
with X (R) and Y (S) amplitudes, characterizing the standard QRPA phonon operator by: 



Y m — T> ym \s m _ At-, ym 

^(pn,J n ) — V Ly pn yV (pn,J 7 ')' 1 (pn,J^) ~ \f u pn I (pn,J*)i 



as: 



R(tt',J-*) — \/'E > tt'R™ t >Jk), S( TT ,j^ — yV TT >S™ T ,^jiry (3.10) 

In the quasiparticle representation, the beta transition density operator can be written 

[CpC n } JM = u p v n A\pn, JM) + u n v p A(pn, JM) 

+ UpU n B^(pn, JM) - v p v n B(pn, JM). (3.11) 

If we restrict our consideration to the ground state to ground state 0^/?/?-decay we end 



up with the following expressions for one-body densities [16,25 



< 



J«m t || [ C ;c n ]j || 0+ > = ^J + l^WxJ^ + vfu®T£^ ) )y/T$ 1 (3.12) 



< 0+ || [c+~c n }j || J*m f > = y/2J + l( v U) v ^X^ >J , ) + u^v^Y^^V^ (3.13) 

Here, the index i (f ) indicates that the quasiparticles and the excited states of the nucleus 
are defined with respect to the initial (final) nuclear ground state |0/" > (|0y >). The overlap 
matrix elements entering Eq. (|2.5|) are explicitly given in Ref. [|39" |. 

The beta transition density from the intermediate states \J' K m > to the first excited + 
state of the daughter nucleus, which is considered to be of two-quadrupole phonon character 
for the nuclei with A=76, 82, 100 and 136, can be written as: 

(0i\[cfd n ]j\J^m>= (0+ PA |-^={r 2 ®r 2 }° {[c + p c n ]j ®Q r $} \0+ PA )V2J + T, (3.14) 

12 



There are two basic approaches to calculate this expression. We shall discuss them in the 
next sections. 

A. The recoupling approach 

The first calculation of the two-neutrino double beta decay (2z//9/?-decay) transition to an 



excited + final state, was presented in Ref. p6[. The formalism proposed was developed in 



the Tamm-Dancoff approximation (TDA). It was claimed that the contribution coming from 
the backgoing graphs are negligible. The dominant contribution is obtained by calculating, 
through a recoupling procedure, the scalar product of two pairs of proton-neutron quasipar- 
ticle creation operators, originating from the beta transition [c+c n ]jM and the phonon Qj M ^ 
operators. 



{A^ipn, J) ® Ai(p>n>, J)}° = - £(-l)-+V+^' ^'[f^ ^f/^1 x 

y 1 (1 + <Vr' I 1 + d nn') ' 

(2J' + 1) 1/2 1 Jn h J \{rf(pp',J')®A\nn',J')Y (3.15) 

[ Jp' Jn' J' J 

Henceforth we shall denote this approach as recoupling method (RCM). 
By using Eq. ( [3.14|) , the beta transition matrix element takes the form 



V U) V U) \ V* 



(Of\{c+c n }j M \.rM,m f >= VlOYsi 1 + V) V2 (1 +^n') 1/2 PP (f r (3-16) 



JSpn 

i Jn Jp J I ^ ^ mf ~p m / -p™/ (/) (/)TF m J' ~Q m f ~Q m f \ 

X \ ( \ U p V n A (p'n',J' r )- K (pp',2+)- K (nn',2+) — V p U n Y (p'n',J' r )'- , (pp',2+)'- , (rm',2+) )' 

jp' jn' 2 J 

Obviously, in the above expression, the full expression of the RPA phonon operator was 
used. In comparing this transition density to 0^ with that one leading to the ground state, 
we find two important differences. First, the dominant contribution in Eq.(3.16) is a product 
of three forward-going amplitudes. This fact implies that the transition amplitude is not 
expected to be very sensitive to the nuclear ground state correlations. Second, the leading 
term in Eq. (3.16) is multiplied by the factor u^f'v^' (i.e., "/3~" like) while the leading term 
of beta ground state transition in Eq. ( |3.13| ) contains the factor v^u)P (i.e., "/3 + " like). 

13 



The drawback of this approach is that the transition density, in Eq. (3.16), contains 
significant unphysical contributions. To clarify this point we transform the second part of 
the r.h.s. of Eq. (3.14), which up to a multiplicative constant should represent the excited 
state 0i~ in the (A,Z+2) nucleus. From the transition operator written in quasiparticle 
representation we keep, for illustration, the operator A ! {pn, J), which is further expressed 
in terms of the pnQRPA bosons. The final result is: 

{A\ P n, J) ® Q^}°Kp Af ) = EWUjraP ® Q$}°K PAf ) 

m' 

^ / o J _i_ 1 Y(£n,J*)\QRPAf)] (3-17) 

The second term in the above equation, is obtained by using the commutator algebra for Q™J 



and its hermitian conjugate operator, and then Eq. (|3.3|) . From Eq. (|3.17| ) it follows that 



within the RCM procedure we have produced a linear combination of a state associated to 
the (A,Z) nucleus and the ground state characterizing the (A,Z+2) nucleus. We note that the 
desired Of excited state of (A,Z+2) nucleus is missing. The RCM does not allow to eliminate 
the admixture of these states, which due to the recoupling procedure are related to the 0+ 
excited state. Moreover, it is worth noting that the component in the (A,Z+2) nucleus 
is proportional to the Y-amplitude, as prescribed by the method presented in the next 
subsection, and not to the forward going amplitude of the proton-neutron dipole phonon 
as suggested by the RCM approach. Thus, the validity of this recoupling procedure is 
questionable in the framework of the QRPA. 

The RCM has been modified by introducing a multiple commutator method (MCM) 

and applied to calculate different lepton-number conserving modes of the double beta decay 

40| . This version of the RCM has been also used for describing the Oz//5/5-decay to excited 



collective + states |23 . 
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B. The boson expansion approach 

The pioneering approach to study the double beta decay to excited states of the final 
nucleus was proposed in Ref. |27],^8|]. It is the so called boson expansion method (BEM). 
Applications of the BEM approach to study the single beta and 2z//3/3-decay to the first 
excited quadrapole state (2+) and the two-quadrapole-phonon states (0^~_ p/l , 2^~_ p/ J of even- 
even isotopes were presented in Refs. p7| , |28| . Recently, the renormalized version of the BEM 



was applied to the transition 82 Se — > 82 Kr |[29| . The new version has the virtue of exploiting 
the complementary features of the BEM and RQRPA methods. As a matter of fact this 
improved version of BEM is adopted in the present paper. 

Within the BEM approach, the operators involved in the r.h.s. of Eq. (|3.11|) are written 



as polynomials of the RPA bosons | 27| , |28 |, so that the mutual commutation relations are 
consistently preserved by the boson mapping. We shall follow this procedure with some 
simplifications, which do not influence the final form of the one-body transition density. 

By exploiting the fact that QjV|0^ Pj4 ) = 0, we introduce a commutator in the expres- 
sion for the one-body transition operator to the 0^ state. In addition, we evaluate this 
commutator by satisfying exact commutation relations. Thus we obtain 

(0t\[c^]jM\rM,m f >= {-l) J - M (^\[[4~ c n]j-M,Q^\Kp A{ ) 

, (f) (f) v m f . (f) (f) v m f s v-^ (°1 I- 8 ( rr 7°0)|0 RPj4 ) , , 

= H < >x (pn,j^ + u P v ny (P lj^) L 77T. — - — -■ ( 3 - 18 ) 

r=p,n J(2j T + 1) 



We omitted the terms A\tt, 00) and A(tt, 00), since in the boson expansion formalism they 
consist of terms comprising products of odd numbers of phonon operators Q j M n , Q ™ M7r , r 2 j u+ 
and T\ M+1 and consequently do not contribute to the above matrix element. 

We proceed by performing the boson expansion of the operator B\tt, 00) with the result 

fit (rr, oo) = Bg(Tr){r 11 ®r&}° 

+ B^{tt){T\ + ® T\ + f + BlKrrnrH ® r^ + } . (3.19) 

The upper indices, accompanying the expansion coefficients, indicate the number of the 
creation and annihilation phonon operators involved in the given terms while the lower 

15 



indices suggest that the phonon operators correspond to the first root of the ppnn-QRPA 
equations. We note that relevant to the problem studied here is the coefficient Bf^rr), 
which can be determined by the following procedure. Commuting Eq. ( |3.19| ) twice with T%+ 
and then taking the expectation value of the result in the boson vacuum, one obtains 

K{rr) = \Y, C 2M 2 -M<0|[r^, K_ M , ^ (tt, 00)]] |0> 

Z M 



— — y 9 • , 1 ( 2^ R{tt' ,2+) S(tt' ,2+) + 2^ R(t't,2+)S(t>t,2+))- (3.20) 

V z Jt T 1 r /( r<r /) t'{t'<t) 

In the above equations, all commutators are exactly evaluated except for the last one for 
which the renormalized quasiboson approximation is used |^7|-^J. 



Now, by a straightforward calculation, one arrives at the final expression for the one-body 
transition density leading to the final excited Of two phonon state: 

(Ot\[cfa]ju\J*M,m, > = (v^uiPxli^+u^^Y^j^ (V pn y 1/2 

€(p,p',n,n'), (3.21) 

where 



£{P,P ,n,n) — Vl0[ — ~( 2^ R(nn',2+)S(nn',2+) + R(n'n,2+)S (n'n,2+)) 

z Jn -T 1 n /( n<n /) 

r, • iil Z^ R(pp',2+)S(pp>,2+) + ^(p'p,2+)^(p'p,2+)) J- (3.22) 

47p "T J- p'(p< P ') 
It is worthwhile to notice that if we replace the factor £(p,p', n, n') with unity and consider 
small ground state correlations, i.e. T> pn ~ 1, we obtain the ground state transition density 
from Eq. ( |3.13|) . We note that the transition density to the ground state is proportional 



to (Upn) 1 / 2 , i.e., it is suppressed by anharmonic effects, while the transition density to the 
excited 0^ state is proportional to l/(T , pn ) 1 / 2 , i.e., it is enhanced by large ground state 
correlations. 



We remark that the BEM expression given in Eq. ( |3.21| ) differs considerably from the 



the RCM expression given in Eq. (3.16). We see that the BEM transition density in Eq. 



( |3.21| ) consists of products of forward and backward going variational amplitudes of the 
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ppnn-RQRPA. It means that the final result exhibits sensitivity to the particle-particle 
interaction of the nuclear Hamiltonian. In addition, the BEM transition amplitude is a 
"/3 + " like amplitude since the leading term is proportional to u^vL This implies a possible 
strong dependence on the ground state correlations. The fact that the RCM and the MCM 
approaches differ considerably from the BEM procedure, when the final state is of multiple 
phonon character, was noticed already in Ref . |28| . The BEM avoids the operators recoupling 



and therefore the problem concerning the unphysical RCM contributions does not appear. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

The formalism described in the previous sections were applied to the transitions 76 Ge — > 
76 Se, 82 Se ^ 82 Kr, W0 Mo ^ 10 ° Ru and 136 Xe ^ 136 Ba. The pn-RQRPA and the ppnn- 



RQRPA calculations have been performed for the same sets of basis-states as in Ref. |T6 
which are identical for protons and neutrons. The single particle energies were obtained 
by using a Coulomb corrected Woods Saxon potential. The realistic interaction employed 
is the Brueckner G-matrix of the Bonn one-boson exchange potential. The truncation 
of the single-particle space requires a renormalization of two-body matrix elements. The 
scaling of the pairing strength in the BCS calculation was adjusted to fit the empirical 



pairing gaps according to Ref. [IlJ. In the RQRPA calculations, the particle-particle and 
particle-hole channels of the G-matrix interaction are renormalized by multiplying them 
with the parameters g pp and g h , which, in principle, should be close to unity. Our adopted 
value for g p h was g p h = 0.8, as in our previous calculations PJl6||. We shall present the 



relevant nuclear matrix elements for g pp = 1. Nevertheless, their sensitivity to g pp within the 
interval 0.80-1.20, which can be regarded as physical, will be discussed. We note that in our 
calculations the pn-RQRPA and ppnn-RQRPA channels are coupled through the equation 
for the renormalization factors D f25f . Our numerical analysis shows that the quadrupole 
QRPA energies of the daughter nucleus are independent of g pp . 

The results of our calculations are summarized in Tables |, [TJ, [TV] and in Fig. [TJ. In Table 
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I the dimensionless nuclear matrix elements of light and heavy neutrino exchange modes of 
0^/9/9-decay of 76 Ge, 82 Se, 100 Mo and 136 Xe, are presented both for transitions to the ground 
and excited states. The displayed Oz//3/3-decay matrix elements to the first excited states 
0^, were obtained within the RCM and BEM approaches. The particular contributions to 
the full matrix elements coming from Fermi, Gamow- Teller and tensor terms in Eqs. Q2.4| ) 
are shown as well. The modifications coming from induced nucleon currents are included in 
the Gamow- Teller and tensor components [16|. We find that the tensor contribution plays 



an important role when the mechanism is mediated by heavy neutrinos and tends to cancel 
the contributions by Fermi and Gamow- Teller transition amplitudes. By glancing at Table | 
we find that the nuclear matrix elements involving the first excited Of state are suppressed 
in comparison with those associated to transitions to the ground state. In the case of BEM 
calculations of M< mv>) the suppression factor is about 2.8, 2.8, 1.8 and 1.5 for A=76, 82, 
100 and 136, respectively. The RCM values are close to the BEM ones for A=76, 82 and 
100 nuclei. One notices an anomaly in the case of the A=136 system, where the RCM 
transition to the excited state is by a factor of 6.7 stronger than that to the ground state. 
It could be connected with the fact that 136 Xe is a closed shell nucleus for neutrons (N=82) 
and therefore the unphysical contributions to this transition in the RCM approach might be 
larger. We note that it is not possible to compare directly nuclear matrix elements Mjy riu> in 
Table |, with those calculated in Ref. |^| for A=76 and 82, since Ref. p3| does not include 



contributions from the induced currents. Nevertheless, we note that the ratio of the nuclear 

matrix elements of the transition to ground and excited states is equal to about 3, in Ref. 

23fl . This value is in good agreement with results of this article, despite the fact that the 



two formalisms differ from each other in many aspects. 

From Table | it follows that the 0z//3/3-decay, mediated by heavy neutrinos to excited 
final states for A = 76, 82 and 100 are weaker than those associated with ground state to 
ground state transitions, by about a factor of 1.3-2.0 ( 8-9 ) in the BEM (RCM). Here, 
the difference between the BEM and RCM predictions is more significant than for the light 
Majorana neutrino exchange mechanism. We note also that for the A=136 system, the RCM 



value of M^ eavy is comparable with the BEM one, i.e., within the RCM, the behavior of this 
nucleus is different from that of the remaining nuclei. 

In Table |TJ, the nuclear matrix elements associated with the trilinear R-parity violating 
mode of Oz//3/3-decay are displayed. Both, the one pion- and two pion-exchange Gamow- 
Teller and tensor contributions to M^, are shown. In Ref. [|3"2"| , |3"3"| , it was shown that there 
is a dominance of the two pion-exchange mode for the Oz//3/3-decay transitions connecting 
the initial and final ground states, due to a larger structure coefficients a 2vr and because of 
strong mutual cancelation of the one-pion exchange Gamow- Teller and tensor contributions. 
We see that the second reason does not hold in the case of transitions to 0^ excited states. 
We have found that the one-pion mode plays a more important role for this transition, giving 
a significant contribution to M x , . By comparing the values of nuclear matrix elements for 
ground and excited state transitions, we see that the second one is reduced by factor of 
2.7-3.4 within the BEM. The RCM values are considerably smaller for A=76, 82 and 100 
systems. A different situation is again found for the Ozz/^-decay of 136 Xe where the RCM 
value is close to the BEM result. 

One purpose of our study is also the sensitivity of results for M<^J>, M^ eavy and M A / m 
to the details of the nuclear model. We have examined the O^/3/5-transition matrix elements 
as function of the renormalization factor for the strength of the particle-particle interaction, 
g pp , considered in the physical interval (see Table |lll|). We see that the BEM values for 



the transition to excited 0^ states exhibit a very similar dependence on g pp as those for 
the transitions to the ground state. On the other hand, the RCM values are insensitive to 
changes of g pp . Thus, our expectations from the previous section, hinging on the forms of 
the BEM and RCM one-body transition densities, have been confirmed. We note that a 
similar behavior has been found also for other nuclear systems. 

As it was already mentioned in the introduction there is additional suppression of the 
0^/9/9-decay to excited + final states coming from the smaller kinematical factor Goi [see Eq. 
Q2.1|)]. The values of Gq\ are given in Table [TV]. One finds that the ratio Goi(0+ s )/G ! oi(0i") 
is about 12, 11, 5.2 and 21 for A=76, 82, 100 and 136 systems, respectively. By this factor 
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are the corresponding half-lifes to excited + final state larger. 

For a given nuclear isotope the characteristics of the 0z//?/9-decay refer both to the nuclear 
matrix element and the kinematical factor. For a chosen isotope, it is worthwhile to introduce 
sensitivity parameters with respect to different lepton number violating parameters. Large 
numerical values of these parameters may define those transitions and isotopes which are the 
most promising candidates for a lepton number violating signal, in the Oz//3/3-decay. These 
parameters are defined as follows |T6|Jl7[| : 



( <m „>(Y) = 10 7 \M <mv> \ ^G 01 year } 
( VN (Y) = lO G \M VN \^G 01 year, 



C x[ll (Y) = 10" \M Kii \^G 01 year. (4.1) 

We listed the sensitivity parameters in Table |y|. By a glance, one finds that within BEM 
the largest sensitivity parameters for the transitions to the Of state are associated with the 
A=100 system followed by the A=82 system . The smallest ones are of the A=76 and A=136 
systems, i.e., these transitions are less favorable for an experimental study. Naturally, there 
are also additional aspects which experimentalists have to take into account in planning a 
search for the Oz//3/3-decay. 

In Table |V[ we present also theoretical half-lifes by assuming < m v >= leV, rj N = 10~ 7 
and Am = 10~ 4 . One finds that in order to get a limit on < m v > of 1 eV, by measuring the 
transition to the excited Of final state, the experimentalists should reach the level of about 
10 28 years for the half-life. By comparing the theoretical values of T®Y 2 for the decay to ground 
state, with the predictions for the transition to the excited state 0^, both yielded by BEM, 
we note that the second ones are larger by about 1-2 orders of magnitude. This situation is 
shown also in Fig. [I]. The question, for experimentalists, is whether the coincidence between 
the deexcitation 7 and the emitted electrons allows to reduce the background in the 0vf3f3- 
decay experiment to a sufficient extent so that the constraints on lepton number violating 
parameters deduced from the transition to the + excited state, can compete with those 
associated with transitions to the ground state. It might be that in the case of the Ovftft- 
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decay of 100 Mo to excited Of state triggered by light or heavy Majorana neutrino exchange 
mechanisms, the suppression of the half-life by a factor of 17 relative to the transition to 
ground state is compensated by diminishing the background events. 

The expected improved experimental upper limits on the 0^/?/?-decay half-life T x J^~ exp , 
imply more stringent limits on lepton number violating parameters < m v >, r] N and A' m . 



By using the sensitivity parameters £'s given in Table [TV], they can be deduced in a straight- 
forward way as follows: 

1(T 5 



< 



m e 



10 24 years 10 



V, 



< 



S<m„> \ i 1//2 ^rj N \ 



10 24 years 

rpOv—exp ' 
i l/2 



/y a2^ 2 / mq \ 4 ( m- g \ 10 r 



10 24 years , . 

100 GeVj V100 GeVj Cx' in \ T°J~ exp ^ ' ^ 

with k = 1.8 |. One finds that in order to push down the upper constraint on < m v > 
below 0.1 eV in the 0^/?/?-decay experiment to excited 0^" final state, one has to measure 
the half-life of 4.02 x 10 28 , 8.96 x 10 27 , 7.59 x 10 26 and 4.77 x 10 28 years for A=76, 82, 
100, 136 isotopes, respectively. Best present limits on this type of decay are on the level 
of 10 21 -10 22 years (see review [^TJ). But we would like to mention that some progress in 



measuring transitions to excited states is expected in the future. For example, experiment 
MAJORANA with 500 kg of 76 Ge plan to have sensitivity of 10 28 years [[y|. 

V. SUMMARY AND CONCLUSIONS 

The 0^/3/3-decay is a sensitive tool to study lepton number violating mechanisms, which 
are associated with the Majorana neutrino mass and R-parity violating supersymmetry. 
Neutrino oscillations indicate that we may be close to the observation of this exotic rare 
process. This stimulates both experimental and theoretical studies, which can be helpful to 
limit the effective Majorana neutrino mass < m u > below the 0.1 eV level. It is an open 
question whether the future 0z//3/3-decay experiments measuring transitions to the excited 
final states can be of comparable sensitivity to different lepton number violating parameters 
as the ground to ground transitions. Experimental studies of transitions to an excited Of 
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final state allows to reduce the background by gamma-electron coincidences. Drawbacks are 
lower Q values and possibly suppressed nuclear matrix elements. The theoretical studies of 
the corresponding nuclear transitions are of great interest. 

We evaluated 0z//3/?-decay nuclear matrix elements for transitions to first excited 0^ fi- 
nal states for A=76, 82, 100 and 136 nuclei. The calculations have been performed within 
two known approaches, the boson expansion method (BEM) and recoupling two pairs of 
quasiparticle operators (RCM). The results of these two types of calculations differ from 
each other considerably especially in the case of exchange of heavy particles for A=76, 82 
and 100 systems. We indicated the drawbacks of the second approach. We also found 
anomalous behaviors of the RCM results for A=136 nuclei. The resulting matrix elements 
are summarized in Table | and Table ||. The suppression of the decay matrix elements to 
Of in comparison with the decay to 0+ s , depends on the isotope and O^/3/3-decay mecha- 
nism. An average suppression factor of about 2-3 is predicted by the BEM. Contrary to the 
RCM results, the BEM ones are significantly depending on the strength of the two-body 
interaction. 

Further, we have calculated sensitivity parameters to different signals of lepton number 
violation associated with transitions to excited final states. We compared them with the 
ground to ground transitions. We have found the largest sensitivity to these parameters 
in the A=100 nuclear system. By comparing with the decay to 0+ we find a suppression 
by a factor of about 4.1 for Majorana neutrino-exchange mechanisms. It means that the 
corresponding theoretical half-life is larger than that associated to the transition to 0+ , by 
factor of 17. In order to reach the sensitivity to a neutrino mass < m v >« 0.1 eV in the 
0^/9/9-decay to an excited + state, it is necessary to measure half-lifes of at least 8 x 10 26 
years (as predicted by BEM). Perhaps, this limit might be expected to be reached by the 
0^/9/9-decay 100 Mo experiment in the near future. 
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TABLES 
TABLE I. Nuclear matrix elements of light and heavy Majorana neutrino exchange modes for 

the O^/3/3-decay in A=76, 82, 100 and 136 nuclei. Both transitions to the ground state (0+ s ) and 

the first Of excited states (which is assumed to be a two-phonon state) of the final nucleus are 

considered. The calculations have been performed within renormalized QRPA with the help of the 

recoupling (RCM) and the boson expansion (BEM) approaches for the evaluation of the one-body 

transition density, (f.s. means "final state" in the table). 





f.s. 


meth. 


Li| 


2;ht neutrino exch. mech. 


Heavy neutrino exch. 
M V F N M& M^ N 


mech. 


A 


M <m v > 


M<™"> 


M< mv> 


M-<m v > 


M v N 


76 


°ts. 




-1.26 


2.18 


-0.190 


2.80 


-37.4 


28.7 


-20.1 


32.6 




of 


RCM 


-0.570 


0.933 


-0.022 


1.28 


-4.00 


1.42 


-0.39 


3.59 




of 


BEM 


-0.371 


0.796 


-0.039 


0.994 


-12.0 


9.62 


-1.06 


16.3 


82 


°ts. 




-1.15 


2.07 


-0.172 


2.64 


-34.4 


25.4 


-17.4 


30.0 




of 


RCM 


-0.617 


0.971 


-0.024 


1.34 


-4.11 


1.40 


-0.014 


4.02 




of 


BEM 


-0.342 


0.762 


-0.033 


0.947 


-11.2 


8.61 


-0.498 


15.2 


100 


<£.. 




-1.28 


2.62 


-0.230 


3.21 


-44.3 


34.2 


-32.9 


29.7 




of 


RCM 


-0.305 


1.059 


0.016 


1.27 


-2.98 


1.21 


0.483 


3.60 




of 


BEM 


-0.397 


1.52 


-0.008 


1.76 


-14.1 


11.1 


-3.90 


16.2 


136 


°ts. 




-0.504 


0.496 


-0.161 


0.66 


-21.7 


16.8 


-16.6 


14.1 




of 


RCM 


-1.66 


3.40 


-0.038 


4.42 


-11.4 


4.37 


0.445 


12.1 




of 


BEM 


-0.205 


0.347 


-0.038 


0.441 


-8.69 


6.98 


-1.99 


10.5 
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TABLE II. The Ou /3 /3-decay nuclear matrix elements associated with trilinear R-parity violating 
mode for A=76, 82, 100 and 136. The same notation is used as in Table I. 





f.s. 


meth. 


ftp 


SUSY mech. 










A 


M& 


Mft 


M l7T 


M& 


M'ft 


M 27r 


My in 


76 


at*. 




1.30 


-1.02 


-24.3 


-1.34 


-0.652 


-601. 


-625. 




of 


RCM 


0.254 


-0.009 


-21.7 


-0.139 


-0.014 


-46.3 


-68.0 




of 


BEM 


0.482 


-0.027 


-40.2 


-0.475 


-0.050 


-158. 


-198. 


82 


o+, 




1.234 


-0.873 


-31.9 


-1.258 


-0.572 


-551. 


-583. 




of 


RCM 


0.253 


0.011 


-23.4 


-0.135 


-0.000 


-40.8 


-64.2 




of 


BEM 


0.462 


0.001 


-40.9 


-0.449 


-0.030 


-144. 


-185. 


100 


°ts. 




1.433 


-1.726 


-25.9 


-1.525 


-1.048 


-775. 


-750. 




of 


RCM 


0.204 


0.018 


-19.6 


-0.102 


0.017 


-25.6 


-45.2 




of 


BEM 


0.532 


-0.208 


-28.6 


-0.509 


-0.129 


-192. 


-221. 


136 


°u 




0.606 


-0.840 


20.7 


-0.742 


-0.543 


-387. 


-367. 




of 


RCM 


0.783 


0.033 


-72.1 


-0.383 


0.021 


-109. 


-181. 




of 


BEM 


0.284 


-0.084 


-17.7 


-0.318 


-0.076 


-119. 


-136. 



TABLE III. The calculated nuclear matrix elements of Oz^/3/3-decay of 76 Ge associated with 
exchange of light and heavy neutrinos and gluinos for different values of g pp within its expected 
physical range in the RQRPA. 







mS ( > 






j^heavy 






My 


in 


9pp 


vts. 


0+BEM 


I RCM 


o+, 


O^BEM 


of 


RCM 


o+, 


0+BEM 


Of .RCM 


0.8 


3.8 


1.34 


1.30 


39. 


19. 






3.7 


-686. 


-225. 


-71. 


1.0 


2.8 


0.99 


1.28 


33. 


16. 






3.6 


-625. 


-198. 


-68. 


1.2 


1.6 


0.58 


1.20 


27. 


13. 






3.3 


-564. 


-163. 


-63. 
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TABLE IV. The sensitivity factors C<m L ,>, Cv , C\ [see Eqs. (4.1)] and calculated O^/3/3-decay 
half-lifes T]_ii for transitions to both ground and excited 0^ states of the final nucleus (A=76, 82, 
100 and 136) by assuming < m u >= leV, rj N = 10~ 7 and Am = 10~ 4 . Goi is the kinematical 
factor. 







7e Ge 


82 Se 


100 Mo 


136 Ie 






0+ -^0+ 
g.s. g.s. 


0v(3(3-dec&y transition 






E t - E f [MeV] 




3.067 


4.027 


4.055 


3.503 


Goi [y- 1 ] 




7.98 x 10~ 15 


3.52 x 10" 14 


5.73 x 10~ 14 


5.92 x 10~ 14 


S<m v > 




2.49 


4.95 


7.69 


1.60 


^N 




2.90 


5.64 


7.10 


3.43 


C\in 




5.57 


10.9 


17.9 


8.92 


T 1/2 (< m u >= 


leV) [yr] 


4.21 x 10 24 


1.07 x 10 24 


4.42 x 10 23 


1.02 x 10 25 


T1/2 (r] N = 10"' 


J ) [yr] 


1.19 x 10 25 


3.14 x 10 24 


1.98 x 10 24 


8.50 x 10 24 


T 1/2 (Am = 10 


~ 4 ) [yr] 


1.04 x 10 25 


2.73 x 10 24 


1.01 x 10 24 


4.07 x 10 24 



®t.s. ~^ 0i Oz^/3/3-decay transition 
Ei - E f [MeV] 1.945 2.539 2.925 1.924 

G i [y' 1 ] 6.58 x lO^ 16 3.25 x 10~ 15 1.11 x 10~ 14 2.81 x 10~ 15 

RCM calculation 

S<m v > 

T-i/2 (< rn v >= leV) [yr] 
T1/2 (r] N = 10- 7 ) [yr] 
T 1/2 (A m = 10- 4 ) [yr] 



•,<m v > 



0.328 


0.764 


1.34 


2.34 


0.092 


0.229 


0.379 


0.641 


0.174 


0.366 


0.476 


0.959 


2.42 x 10 26 


4.47 x 10 25 


1.46 x 10 25 


4.76 x 10 24 


1.18 x 10 28 


1.90 x 10 27 


6.95 x 10 26 


2.43 x 10 26 


1.06 x 10 28 


2.42 x 10 27 


1.43 x 10 27 


3.52 x 10 26 


BEM calculation 






0.255 


0.540 


1.85 


0.234 
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C'/, 



c 



An 



T1/2 (< m u >= leV) [yr] 
T1/2 (r] N = W- 7 ) [yr] 



r i/2 (^111 



10" 



[yr] 



0.418 

0.508 
4.02 x 10 26 
5.72 x 10 26 
1.26 x 10 27 



0.866 

1.055 
8.96 x 10 25 
1.33 x 10 26 
2.91 x 10 26 



1.71 
2.33 

7.59 x 10 24 
3.43 x 10 25 
5.97 x 10 25 



0.557 

0.721 
4.77 x 10 26 
3.23 x 10 26 
6.23 x 10 26 
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FIG. 1. Calculated half-lifes of the Ov (3 /3-decay of 76 Ge, 76 S"e, 100 Mo and 136 Xe for transitions 
to the ground 0^ s and 0J~ excited states of the final nuclei assuming < m v > ee = leV (a), 77^ = 10~ 7 
(b) and A' m = 10 (c). The black bars correspond to results describing ground state to ground 
state transitions as well as the results obtained for the transitions to the Of excited state within 
the boson expansion method (BEM). The open bars denote results obtained for transitions to the 
0± state via the recoupling method (RCM). 
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